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We study an undoped t-J model with electron-phonon interaction using the self-consistent Born 
approximation (SCBA). By neglecting vertex corrections, the SCBA solves a boson-holon model, 
where a holon couples to phonons and magnons. Comparison with exact diagonalization results for 
the t-J model suggests that the SCBA describes the electron-phonon interaction fairly accurately 
over a substantial range of J/t values. Exact diagonalization of the boson-holon model shows that 
the deviations are mainly due to the neglect of vertex corrections for small J/t and due to the 
replacement of the t-J model by the boson-holon model for large J/t. For typical values of J/t, 
the electron-phonon part Sep of the electron self-energy has comparable contributions from the 
second order diagram in the electron-phonon interaction and a phonon induced change of magnon 
diagrams. A very simple approximation to Eep gives a rather accurate effective mass. Using this 
approximation, we study the factors influencing the electron-phonon interaction. Typically, we flnd 
that the magnons nominally have a stronger coupling to the holon than the phonons. The phonons, 
nevertheless, drive the formation of small polarons (self-localization) due to important differences 
between the character of the phonon and magnon couplings. 



I. INTRODUCTION 

There have recently been several experimental indica- 
tions that the electron-phonon interaction plays a sub- 
stantial role for properties of cuprates, for instance in 
photoemissioniSii and neutron scattering workj^ The 
strong effects of the Coulomb interaction in the cuprates 
is often taken into account by using the t-J model. ^ In- 
cluding phonons in the t-J model, it was concluded that 
the Coulomb interaction can enhance the effects of the 
electron-phonon interaction for undoped cuprates. ^'^ It 
is then interesting to study this aspect further. A simple 
method for treating the undoped t-J model is the self- 
consistent Born approximation (SCBA) (SiSiiSiiiiiS which 
can also be applied to the t-J model with phononsiSii^ 
This method assumes a quantum Neel ground-state for 
the undoped system. The excitations of the spin sys- 
tem are described by antiferromagnetic magnons. A hole 
created in, e.g., photoemission is assumed to interact 
with magnons and phonons, which are both treated as 
bosons. This model is here referred to as the boson- 
holon model. The electron self-energy of this model is 
expressed in terms of the simplest diagrams, including a 
boson (magnon or phonon) Green's function and a self- 
consistent electron Green's function. Vertex corrections 
are neglected. 

Here we extend the SCBA study of Ramsak et al.,^ fo- 
cusing on the limit of weak electron-phonon coupling for 
0.2 < J/t < 2. For strong electron-phonon coupling the 
SCBA is known to break downi We first compare results 
for the quasiparticle weight and energy using the SCBA 
and exact diagonalization of the t-J model. The results 
suggest that the SCBA describes the electron-phonon in- 
teraction reasonably well for a substantial range of J/t 
values, but that it is less accurate than for the t-J model 
without phonons. To trace the sources of errors in the 
SCBA, we use exact diagonalization for the boson-holon 
model. Comparison between these results and results 
from the SCBA shows that for small J/t(^ 0.2) errors 



in the SCBA are mainly due to the neglect of vertex 
corrections in the SCBA, while for large J/t the main 
source of errors is the replacement of the t-J model by the 
boson-holon model. We then study the electron-phonon 
contributions to the electron self-energy. There is a con- 
tribution from the diagram containing one phonon and 
one electron Green's function, which for noninteracting 
electrons is the leading contribution. Here there is a com- 
parable contribution from diagrams containing magnons 
and one electron Green's function due to the change of 
the self-consistent electron Green's function induced by 
the electron-phonon interaction. We also study the effec- 
tive mass. By slightly modifying a previous approaches 
we obtain a very simple formula for the effective mass, 
which agrees rather well with exact results within the 
SCBA. This formula is used to illustrate the factors in- 
fluencing the strength of the electron-phonon coupling in 
the undoped t-J model. We discuss the important differ- 
ence between the coupling to magnons and phonons in 
terms of strength and effects of vertex corrections. We 
comment on the implications for formation of small po- 
larons (self- localization). 



The boson-holon model and the SCBA are described 
in Sec. The SCBA results are compared with exact 
diagonalization results for the t-J model in Sec. IIIII In 
Sec. IIVI we compare with exact diagonalization for the 
holon-boson model to determine the sources of errors in 
the SCBA. The contributions to the electron-phonon part 
of the electron self-energy and to the effective mass are 
discussed in Sec.lVland Sec. I VII respectively. In Sec. I VIII 
we compare the coupling to magnons and phonons and 
discuss polaron formation. 
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II. MODEL AND METHOD 

The t-J modeiS, is given by the Hamiltonian 

Ht-j = J ' — Y 

<ij> 
<i,j>a 



(1) 



where cj^ creates a hole on site i if this site previously 
had no hole. The Zhang-Rice singlets are represented 
by empty sites. Here t is a hopping integral, J is the 
exchange interaction, is the spin on site and 
is the occupation of site i. We use the electron-phonon 
interaction 



(2) 



where N is the number of sites and 6]^ creates a phonon 
with the wave vector q. We assume an on-site coupling 
with the strength ^q. The coupling to hopping integrals 
and to the spin-spin interaction are neglected, since these 
couplings have been found to be weak^ In the following 
we in addition assume a Holstein type of couphng, i.e., 
gq = g is q-independent. 

Following previous worl&SiSiiSiiiii^ for the undoped 
system, the Hamiltonian Ht-j + H^p is approximately 
rewritten in terms of a boson-holon model, where spin- 
less holons interact with phonons and antiferromagnetic 
magnons, treated as bosons, 

kq 



(3) 



where h^^ and aj^ create spinless holons and antiferro- 
magnetic magnons, respectively. The fermion-magnon 
coupling is given by 



Mkq = VSt 



7k-q 



-l-7kSgn(7q)A/ J/q^-l 



(4) 



where 7q = {cosq^ + cosqy)/2 and i/q = (1 — 7q)^^^. The 
magnon frequency is given by ujq = 2Jvq and the phonon 
frequency by Uph- 

The Hamiltonian H is treated in the self-consistent 
Born approximation. The electron self-energy is then 
given by6i8i9ii°iiiii2 

E(k,c.) =l^[Af2qG(k-q,c.-c.q) (5) 



-S'qG'(k-q,tj-tjph)], 



where G'(k, uj) is the holon Green's function, 

1 



G(k,c.) 



CO — S(k, Lo) 



(6) 



Putting ^q = 0, we obtain the corresponding quan- 
tities without electron-phonon coupling, Go(k, w) and 
Eo(k, w). We also introduce the electron-phonon part 
of the electron self-energy 



I]op(k,w) = i;(k,a;) - i;o(k,w), 
and split Eop(k,w) in two contributions 

I]2nd(k, c.) = 1 ^ 4G(k - q, - C^ph) 



(7) 



(8) 



and 



^^piKcu) =l^[Af^JG(k-q,c.-c.q) (9) 
q 

-Go(k-q,w-tjq)]. 

Here I]^p'^(k, ui) corresponds to the second order diagram 
in the electron-phonon coupling. For noninteracting elec- 
trons this is the leading contribution in to S. For the 
interacting system there is a second contribution of the 
same order in g, E^(k, lj). This is due to the change 
of the contribution from the diagram describing the cou- 
pling to magnons caused by to the change of the Green's 
function when the electron-phonon coupling is turned on. 
We also introduce^ 



S,Cp°nk,c.) = - ^5^G[r°'^(k - q,c. - ..ph), (10) 



N 



where GQ°'^(k — q, w — Wph) only includes the coherent 
part of the Green's function 



G^°^k,c.) = 



UJ - eo(k) ' 



(11) 



Here ^o(k) and £o(k) are the quasiparticle strength and 
energy, respectively, in a system where g = 0. Since we 
consider the limit of weak electron-phonon coupling be- 
low, we have neglected the effect of the electron-phonon 
coupling on Z(k) and £(k) in Eqs. (|10l The quasi- 
particle energy is determined by the Dyson equation 



£o(k) = Eo(k,eo(k)). 



(12) 



The shift of the quasiparticle energy due to the electron- 
phonon interaction is then 

Ae(k) = e(k) - £o(k) ~ Zo(k)S]ep(k, £o(k)), (13) 

where 



^o(k) 



1 



aso(k,c^) 



w=eo(k) 



(14) 
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FIG. 1: Az/Ao [Eqs. JTSHHJ] for a;ph/t = 0.1 as a function of 
J/t for a 4 X 4 cluster according to exact diagonalization of 
the t-J (full line) and the boson-holon (dashed line) models, 
the SCBA (dotted line) and SCBA together with the lowest 
order vertex corrections (dash-dotted line) in the limit of a 
small coupling g. 



III. COMPARISON WITH EXACT 
DIAGONALIZATION 

There have been extensive comparisons of results from 
the SCBA and exact diagonalization for small clusters 
for the case of no electron-phonon intcractionjSiSiiSiiiiiS 
and the two methods have been found to agree rather 
well. Here we therefore focus on the changes due to the 
electron-phonon interaction. We define a dimcnsionlcss 
electron-phonon interaction Xz from 



Zo(V2,7r/2) 
Z(7r/2,V2) 



l = Ap 



(15) 



where Z(7r/2, 7r/2) is the quasiparticle weight for k = 
(7r/2, 7r/2) including the electron-phonon coupling. We 
study a 4 X 4 lattice with periodic boundary conditions 
in the limit of weak electron-phonon coupling, for which 
exact diagonalization can easily be performed. Following 
earlier work,^-^ we define Aq as the corresponding quan- 
tity for a single electron at the bottom of the band of 
a two-dimensional Holstein model with nearest neighbor 
hopping. Assuming Uph/t <C 1 and using a quadratic 
expansion of the band, we obtain 



Ao 



ph 



(16) 



We emphasize that by considering the bottom of the band 
the resulting Aq is particularly smalL"'^- At larger filling 
the corresponding Aq is larger and the resulting enhance- 
ment of Xz is smaller. It is necessary to pay some extra 



attention to the q = coupling, go- By using Eqs. (|10l 
111! I13|l . we find that in the SCBA this component gives 
a contribution 



1 

N 



90 



(17) 



to Xz- The q = component just couples to the total 
number of electrons. Since the Green's function describes 
the addition or removal of an electron, we can alterna- 
tively calculate the exact spectrum for one electron cou- 
pling to the q = component and then convolute this 
spectrum with the spectrum resulting from the coupling 
to the q 7^ components. The exact q = contribution 
to Xz is then 



N Vwph/ 



(18) 



This diflFers from Eq. l(T7|) by a factor of Zo(k)^, which is 
typically a very large difference. Although the compar- 
ison with exact diagonalization can only be done for a 
small cluster, we are primarily interested in infinite sys- 
tems where the q = component plays no role. In dis- 
cussing Xz and Ae below we therefore exclude the q = 
coupling. 

FigureQlshows exact results (full line) and results from 
the SCBA (dotted line) for iVph/t = 0.1. The dashed 
and dash-dotted curve are discussed in Sec. IIVI The 
results agree qualitatively. However, quantitatively the 
agreement is not as good as found for the model without 
electron-phonon interaction. This is, in particular, the 
case for small J/t, which are values usually assigned to 
the high- Tc cuprates, and for large values of J/t. For 
instance, Xz/Xq is about 3.3 and 5.2 according to the ex- 
act calculation and the SCBA, respectively, for J/t = 0.3. 
Figure 13 shows the energy shift Ae(7r/2, 7r/2) [Eq. ifT^ ] 
due to the electron-phonon interaction for a 4 x 4 clus- 
ter. As in the case of Xz, the SCBA deviates appreciably 
from the exact results for small and large values of J/t. 
The reason for these deviations are discussed in Sec. IIVI 
The agreement with exact results is still good enough 
to suggest that we can use the SCBA for a qualitative 
discussion of properties of the t-J model with phonons. 



IV. ACCURACY OF APPROXIMATIONS 
BEHIND THE SCBA 

In view of the results in Fig. ^ and Fig. |21 it is in- 
teresting to ask for the sources of errors in the SCBA. 
We distinguish between two classes of errors: i) the re- 
placement of the t-J model in Eq. ^ by the boson-holon 
model in Eq. and ii) the neglect of vertex corrections 
when solving this model. 

We consider the first class of errors by solving the 
boson-holon model [Eq. (Q] using exact diagonalization 
in the limit of weak electron-phonon coupling. The 
Hilbert space can then be limited by only considering 
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FIG. 2: Energy shift Ae for k = (7r/2, 7r/2) [Eq. lO] for a 4 x 
4 cluster with tJph/i = 0.1 according to exact diagonalization 
of the t-J (full line) and the boson-holon (dashed line) models, 
the SCBA (dotted line) and SCBA together with the lowest 
order vertex corrections (dash-dotted line) in the limit of a 
small coupling g. 



a) b) 

Magnoii. Magnon Magnon, ^, , Phonon 



FIG. 3: Self-energy diagrams including (a) two crossing 
magnon lines and (b) crossing magnon and phonon lines. 
The magnon, phonon and self-consistent electron propaga- 
tors are represented by dashed, dash-dotted and full lines, 
respectively. 



states with at most one phonon excited. Due to the 
strong holon-magnon coupling, however, it is necessary 
to consider states with many excited magnons. Figures^] 
andUl compare these results with results from diagonal- 
izing the t-J model. Except for large J/t, the agreement 
is very good. This shows that for small and intermediate 
values of J/t, the errors in the SCBA are mainly due to 
vertex corrections, while for large J/t the replacement of 
the t-J model by the boson-holon model leads to appre- 
ciable errors. 

We next consider vertex corrections. Figure|3Ji shows a 
second order diagram which is not included in the SCBA, 
due to the crossing magnon lines and resulting vertex 
correction. It was shown by Liu and Manousakis^- that 
this diagram and many other diagrams neglected in the 
SCBA are actually zero due to the symmetry of Mkq. 
This makes it understandable why the neglect of vertex 
corrections is found to be a rather good approximation 
for magnons. 



In the limit of weak electron-phonon coupling discussed 
in this paper, diagrams with crossing phonon lines do 
not contribute to lowest order in g^. Diagrams involv- 
ing magnon line(s) crossing one phonon line, however, 
do contribute to this order in g^. Figure shows the 
lowest order diagram of this type. In contrast to the 
pure magnon diagram in Fig. the diagram in Fig. 
is in general not zero. We have included this diagram 
and an equivalent diagram in the calculations, using self- 
consistent propagators for all electron lines in the calcu- 
lation of Xz and Ae. The results are shown in Figs. 
and 121 by the dash-dotted lines. For small J/t the cor- 
rection to the SCBA (dashed curve) is large (almost a 
factor of two) and it goes in the correct direction com- 
pared with the exact result for boson-holon (dashed line) 
and t-J models (full curve). Some higher order diagrams 
are not small, although the sum of all higher order dia- 
gram apparently almost cancel. The rapid convergence 
for small J/ 1 suggested by Fig. and Fig. El is therefore 
somewhat misleading. We note that these lowest order 
vertex corrections have a substantially smaller effect on 
the self-energy for larger clusters. In these cases, how- 
ever, it is not possible to perform exact diagonalization 
calculations, and it is therefore not clear if the SCBA 
becomes more accurate for large clusters. 

We observe that diagrams of the type in Fig. were 
neglected in the calculation of the criterion for polaron 
formationpi If the results in Fig. and Fig. (21 can be 
extrapolated to large clusters and strong coupling, they 
suggest that the earlier criterion^ for polaron formation 
may have underestimated the critical A. 

We have elsewhere studied vertex corrections to the 
electron-phonon interaction for the half-filled Hubbard 
model in the large U limit, which is closely related 
to the undoped t-J model. Neglecting vertex correc- 
tions and considering weak electron-phonon coupling, we 
found that a sum rule for the electron-phonon part of 
the imaginary part of the electron self-energy is strongly 
violatedfi^ in apparent contradiction to the fairly good 
results found in the SCBA above. The violation of the 
sum rule in the large-?/ Hubbard model can be traced to 
the fact that the weight of the spectral function only in- 
tegrates to one half over the photoemission energy range. 
This problem is avoided in the SCBA by using the spin- 
less holon Green function, for which the spectral func- 
tion integrates to unity (over the photoemission energy 
range) . 



V. QUASIPARTICLE ENERGY 

In the remainder of the paper we focus on the SCBA 
and first analyze the quasiparticle energies. We have 
performed calculations for 96 x 96 lattices, using t = 1, 
J/t = 0.3, g/t = 0.05 and Wph/i = 0.1. The self-energy 
was broadened by adding a small imaginary part S/t — 
0.005 — 0.01 to the energy. Figure 01 compares results 
for Ae(k) (full line), determined from the Dyson equa- 
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FIG. 4: Zo(k)Ecp(k,£o(k) as a function of k along the 
(0, 0) - (7r/2, 7r/2), (7r/2, 7r/2) - (0, tt) and (0, tt) - (0, 0) direc- 
tions. Result are shown for [Eq. ^ (dashed hne), Sep"* 
[Eq. ^] (dotted line) and Ecp [Eq. 10] (dash-dotted hne) as 
well as for the exact Ae(k) (full line). The inset shows the 
dispersion of eo(k). The parameters are t = 1, J/t = 0.3, 
g/t = 0.05 and ojph/f = 0.1. 

tion using the full self-energy, with Zo(k)Sop(k, £o(k)) 
[Eq. using three approximations for the self-energy. 
The figure illustrates that S^p (dashed line) is a rather 
good approximation to Epp'^ (dotted line), i.e., the inco- 
herent part of Go included in Sgp'* does not contribute 
much to the self-energy. It is interesting, however, that 
E^p*^ (dotted line) is a rather poor approximation to Eop 
(dash-dotted line) , and it only contributes about half the 
magnitude for J/t = 0.3. Both S^^'i and Efp are of the 
order g^. For large values of J/t this difference is smaller. 
For noninteracting electrons Sop is the only contribution 
of this order, and the interest has therefore often focused 
on this contribution. Figure 01 hows that this is not a 
good approximation for the present model and J/t = 0.3. 
The full line and the dash-dotted line differ slightly since 
the solution in Eq. (|13|l of the Dyson equation is only 
approximate for a finite g. 

To better understand the results for Sop'^; notice 
that for a q-independent = g, Im E^p'^ takes a very 
simple form 

ImE^f (k, u;) = ng'Aiu; - c.ph), (19) 

where A{uj) ^ Y^j^lmGo{k,uj - iO+)/{NTT) is the k- 
averaged spectral function. Figure [S] shows A{uj). Since 
we used cuph/t = 0.1 in Fig. 0] the onset of Im Sep 
has been shifted by OAt above the bottom of the band. 
States below this onset are then shifted strongly down- 
wards, while states above the onset are shifted less or 
are even shifted upwards. From the inset of Fig. 01 we 
can see that states around {tt/2, 7r/2) and along the line 
(7r/2,7r/2) — (tt, 0) are below the onset and are shifted 
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FIG. 5: Spectral function A{uj) with (full line) and without 
(dashed line) electron-phonon coupling. The parameters are 
the same as in Fig. 21 except g — O.lt. A Lorentzian broad- 
ening with FWHM 0.2t was used. 



strongly downwards, in particular states which are just 
below the onset, while Sgp'^ becomes positive for k- 
vectors along the lines (0, 0)-(7r/2, 7r/2) and (tt, 0)-(0, 0) 
close to (0,0). 

In a similar way, we can understand S^p [Eq. ®], al- 
though in this case the coupling Mkq and the energy 
ujq have strong k- and q-dependencies. S^ is due to 
the coupling to the changes of A(k,uj) caused by the 
electron-phonon coupling. Figure |31 shows the k-average 
A{uj) with (full line) and without (dashed line) electron- 
phonon coupling. The parameters are the same as for 
Fig. 21 except that g/t ~ OA to enhance the effect of the 
electron-phonon coupling. This difference in A(uj) is pos- 
itive and particularly large at « —2.2. Since lu^ can 
be fairly large, ranging from zero to 2J — 0.6t, Im S^ 
is shifted substantially upwards in frequency. As a result 
Re S^ is negative for the whole quasiparticle band. 

VI. EFFECTIVE MASS 

We next consider the effective mass, essentially follow- 
ing Ramsak et al.^ Since S^ is relatively k-independent, 
we neglect it and we only consider S^"*^. For simplicity, 
we assume that Zo(k) = Zo(7r/2, 7r/2) is k-independent. 
We furthermore approximate the quasiparticle disper- 
sion by assuming that it can be expanded quadratically 
around the four minima (±7r/2, ±7r/2). Two masses are 
introduced, m\\ and m±^, which describe the dispersion 
parallel and perpendicular to the (0, 0) — (tt, tt) direc- 
tion, respectively. The summation over the Brillouin 
zone in Eq. H10|) is replaced by an integration over all 
of q-space, assuming that contributions far away from 
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(±7r/2, ±7r/2) are small because of the large energy de- 
nominator in Eqs. UlUI By using the solution of the 
Dyson equation Eq. H13|l. we then obtain 

^ ' \2ti) J £o(k)-£o(k-q)-cjph 

(20) 

where the factor four is due to the presence of four equiva- 
lent minima (±7r/2, ±7r/2). One factor of Zq comes from 
the Green's function and a second factor from solving 
the Dyson equation. Defining the effective mass along 
the parallel direction as \/m* = (i^e(k)/(ifc| , we obtain 

where the second equality defines the electron-phonon 
coupHng Am. Focusing on large J/i, Ramsak et al.^ ob- 
tained the same result except for a factor Zq resulting 
from the Dyson equation Without this factor the 

rather good agreement with the exact result in Fig. 
would be lost for small J /t. \,n is compared with the 
corresponding quantity for the Holstein model(& which is 
identical to the Ao defined via Z in Eq. (|16l) . 

Figure ini shows results for Am/Ao as a function of J/i, 
using the second derivative of the exact £(k) (full line) 
and of e(k) obtained from (dotted line) as well as 

Eq. (dashed fine). Am/Ao is different from A^/Aq in 
Fig. n The main reason for this difference is that Fig. 
shows result for a 4 x 4 cluster while Fig. El shows re- 
sult for a large cluster (96 x 96 or 128 x 128), but the 
two quantities are somewhat different also for identical 
clusters. The results based on (dotted line) agree 

rather well with the exacts result for intermediate and 
large values of J/t, while they are too small for small 
values of J/t. The deviation for small J/t is primarily 
due to the neglect of E^. This term is particularly im- 
portant for small J/t, since the magnon energy entering 
in Eq. Q is proportional to J. For large J/t, the devi- 
ation is mainly due to the neglect of the incoherent part 
of G. Eq. 1)2 l|l (dashed curve), which is an approxima- 
tion to the dotted curve, gives a larger A„i and it agrees 
better with the exact result. For small J/t, the increase 
is primarily due to neglect of k-dependence of Z^i/k) in 
the Dyson equation when deriving Eq. H21|) . For large 
large J/t the increase is due to several small errors in 
the approximations. The coupling g /t = Q.Ws somewhat 
too large to give the weak-coupling limit, in particular 
for large J/t. 

We are now in the position to interpret the enhance- 
ment of Am compared with Aq obtained at the bottom 
of the band for a Holstein model. Equation H21I) con- 
tains a factor Zq which tends to reduce the coupling due 
to the transfer of spectral weight far away from the en- 
ergies studied. On the other hand, the factor ^m\f^m_\_ 
describes how the energy denominator is reduced by the 
large effective masses, bringing spectral weight closer to 
the relevant energies;^ This is an effect of correlation and 
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FIG. 6: A,„/Ao as a function of J/t calculated from the quasi- 
particle energy e(k) (full line), Eq. I21II (dashed line) and 
using the self-energy Ejpp'^ [Eq. llOt l. The parameters are the 
same as in Fig. 2] except that g = O.lt and J/t is varied. 

antiferromagnetism, and it is important for the enhance- 
ment of the electron-phonon interaction in this model. In 
addition, there is a factor of four resulting from the pres- 
ence of the four equivalent minima (±7r/2, ±7r/2). For 
J/t = 0.2 we find that Z^ — 0.05 and ^m\\m± — lOmo, 
where mg = l/(2i) is the mass at the bottom of the 
band in the Holstein model. In this case the factor 
four from the equivalent minima is crucial, since the 
electron-phonon interaction would otherwise have been 
suppressed in the t-J model, while now it is enhanced by 
a factor of 1.8. For J/t = 2, we obtain Zq = 0.57 and 
^m\\m±_ ~ 6.5mo, giving the enhancement 16 (15 ac- 
cording to Eq. (|21|l ). In this case, the large mass plays a 
crucial role for the enhancement of the electron-phonon 
interaction. We notice, however, that the comparison 
here has been done with Aq calculated for a single electron 
at the bottom of the band of a Holstein model. Had the 
comparison been made with a half-filled Holstein model, 
the result would have been a smaller enhancement or no 
enhancement at all>i^ By starting from the three-band 
model and studying the half-breathing phonon, however, 
it is found that the coupling constants (7q are enhanced 
by correlation effects 

VII. COMPARISON OF COUPLING TO 
MAGNONS AND PHONONS 

We define an average dimensionless coupling constant 
for the magnons 

Am = ^ 2^ AMk - ^ 2^ ^ - ^ , (22) 
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where AM(7r/2,7r/2) = 0.65^/ J. For La2Cu04 the cor- 
responding quantity due to phonons is A = 1.2.^^ For 
a typical value J/t = 0.3fi^ the coupUng to magnons, 
Am = 1-67 and AM(7r/2.7r/2) = 2.2, is stronger than the 
coupling to phonons. It might then seem that the cou- 
pling to magnons is more important for the experimen- 
tally observed^ polaron formation in undoped cuprates. 
This is, however, misleading. The value of A needed for 
formation of small polarons is reduced with the boson 
frequencyiSSiSiiSS This somewhat favors phonons, since 
they typically have lower frequencies than the magnon 
frequencies of the order of J. To see the main difference, 
however, it is necessary to consider vertex corrections. 

To describe the formation of small polarons due to 
phonons, it is crucial to go beyond the SCBA, since 
vertex corrections including phonon propagators become 
very important in the strong-coupling limit i Actually, if 
these vertex corrections are neglected, polaron formation 
is not properly obtained.^ On the other hand, it has been 
argued that vertex corrections including magnon propa- 
gators are not very important in the t-J modeliSiSiiiiiS As 
discussed in Sec. lIVI the lowest order vertex correction in 
Fig. is identically zero due to the symmetry of Mkq, 
and classes of higher order vertex corrections are also 
zero.^^ If we assume a Holstein type of electron-phonon 
coupling, however, there are no similar arguments for 
diagrams with crossing phonon lines being zero. This 
explains why the holon-phonon interaction, but not the 
holon-magnon interaction, leads to polaron formation. 



VIII. SUMMARY 



To summarize, we have studied the self-consistent Born 
approximation (SCBA) in the limit of weak electron- 
phonon coupling. While the SCBA has been shown to 
be quite accurate for a pure t-J model, we find that it is 
less accurate when the electron-phonon interaction is in- 
cluded. To study the reason for this, we performed exact 
diagonalization calculations for the boson-holon model. 
Comparing the results with the SCBA results, we find 
that the main errors of the SCBA are due to the neglect 
of vertex corrections for small J/t and due to the the in- 
troduction of the boson-holon model itself for large J/t. 
Studying the electron-phonon part of the self-energy, we 
find that in addition to the second order term S^p'* known 
from the theory of noninteracting electrons, there is a 
second term of the same order, S^. For J/t w 0.3 this 
term makes a similar contribution to Sop as S^p^^. We 
have shown that a very simple derivation of the effec- 
tive mass gives a rather accurate result, illustrating the 
factors enhancing and suppressing the electron-phonon 
coupling. The coupling to magnons can be considered 
stronger than the coupling to phonons for realistic pa- 
rameters. Nevertheless, the phonons drive the formation 
of small polarons for undoped cuprates, due to the dif- 
ference between phonons and magnons, in particular the 
different importance of vertex corrections. 
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